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Abstract 


For PEM fuel cells supplied with air, pressure and flow control is a key requirement for an efficient and dynamic operation because fuel cells 
are in risk of starvation when the partial pressure of oxygen at the cathode falls below a critical level. To avoid oxygen starvation and, at the same 
time, to allow for a dynamic operation of the fuel cell system, both excess ratio of oxygen and cathode pressure need to be adjusted rapidly. 

In this paper a model-based control structure is proposed that comprises a multivariable control of the cathode pressure pc and the excess ratio 
of oxygen Ao, using the mass flow controller (MFC) and the outlet throttle as actuators. For both controlled variables appropriate desired values 
can be specified independently. As the derived dynamic model of the air supply system is nonlinear, advanced control techniques using differential 


flatness are applied. 


Since the partial pressure of oxygen in the cathode flow field is a state variable that is hardly accessible by measurement, a tracking observer is 
employed to estimate this pressure using the measured air pressure at the outlet throttle. Thereby, the observer enables the control system to detect 
both temporary shortages of oxygen and peak pressures at the cathode and to counteract accordingly. 


© 2007 Elsevier B.V. All rights reserved. 


Keywords: Starvation prevention; Pressure control; Air flow control; Pressure observation; Fuel cell dynamics 


1. Introduction 


Fuel cells (FC) are energy converters that are not subject to 
the Carnot cycle and can, therefore, provide electrical power and 
heat efficiently, even in the part-load range. In case of hydrogen 
as fuel no direct emissions occur. 

In automotive and several niche transportation applications 
FCs act as propulsion engine or as auxiliary power unit. In sta- 
tionary applications fuel cells are utilised as combined heat and 
power plants or as uninterruptible power supplies. 

The prevalent type of fuel cell is the polymer electrolyte mem- 
brane fuel cell (PEMFC), which is characterised by a high power 
density and a low operating range with temperatures from 0 to 
90°C. In both applications mentioned, a dynamic and efficient 
performance is required due to varying operating conditions and 
an operation with alternating loads. The power response of the 
fuel cell system depends mainly on the interaction of the fuel cell 
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stack itself and the necessary peripheral components like air and 
fuel supply, flow and pressure control, heat and water manage- 
ment. Hence, the control objectives are to maintain efficiency 
and to avoid degradation by guaranteeing appropriate tempera- 
ture, humidity and partial pressures of the reactants across the 
electrode. 

Flow and pressure control in particular is a key requirement 
for an efficient and dynamic performance of the PEMFC [6]. 
Fuel cells are in risk of starvation during high current demand 
(see Fig. 2) and fast load changes (see Fig. 1) [12]. Starva- 
tion comes along with a drop in partial pressure of oxygen and 
occurs when more gas molecules react than are fed to the active 
surface. 

The mass flow mair at the inlet of the cathode is usually 
provided by a feedforward control proportional to the actual 
current, according to Faraday’s Law 


MQp,react 


N Mo, lfc 
= ro 


2 (1) 
YOz 7 YO, 2 F 


Mairin = 407 


with the mass fraction of oxygen in air yo, and the excess ratio 
AO, = M0,,in/MO,,react- AS can be seen in Fig. 1, this feedfor- 
ward control of the mass flow leads to a temporary undersupply 
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Fig. 1. Excess ratio of oxygen and cathode pressure at current step. 


with ào, < 1. Due to the time constant of the supply system 
oxygen cannot be replenished fast enough. 

Diffusion is a limiting factor in the case of high current den- 
sity. Reactant concentration drops and concentration overvoltage 
increases significantly due to both high reactant consumption 
and head loss at high flow rate. The limiting current with every 
incoming molecule reacting immediately is referred to as Jim. 
The limiting current leads to a steep decline in voltage in the IV 
characteristics at high currents (Fig. 2). 

The J-V characteristics can be described theoretically (2) by 
the equilibrium voltage via the Nernst equation and additional 
loss terms in terms of overvoltages nj; [1]. The fuel cell voltage 
Urc is a function of the fuel cell current Ipc, the stack temper- 
ature Trc, the partial pressures po, and py, of the reactants at 
the electrodes, the porosity € of the gas diffusion layer and the 
conductivity « of the membrane 


Urc(Irc, x) = Ueg — Nact — Neone — Nohm; (2) 


with x = [Tfc, Po, Pro, €, K|". 

Related to the partial pressures of the reactants, two effects 
are important. First, decreasing partial pressures lead to a lower 
voltage level all over the J-V characteristics due to an increased 
concentration overvoltage. Second, the limiting current Jim 
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Fig. 2. I-V characteristics with limiting current. 


depends on the stack temperature, the reactant partial pressure 
and the porosity of the gas diffusion layer 


Tim = f(p, T, ©). (3) 


If the partial pressures of the reactants drops, [jim decreases and 
the J-V characteristics is deformed towards smaller currents. 
This reduces the available operating range. A drop in partial 
pressures at high currents puts a crucial risk to the cell. If the 
limiting current falls below the actual fuel cell current, the cell 
polarity gets reversed. This effect harms the cell by damaging 
the catalyser. 

In order to prevent starvation, feedforward control of the air 
mass flow and pressure control at the outlet of the system are 
insufficient. As will be shown in the following, multivariable 
feedback control of the cathode pressure pc and the excess ratio 
of oxygen Ao, using the mass flow controller (MFC) and the out- 
let throttle as actuators, complies with the requirements. Since 
the flow dynamics is nonlinear, as described in Section 3, nonlin- 
ear model-based pressure control is proposed in Section 5. Here, 
a first problem arises: the partial pressure of oxygen at the elec- 
trode is an unknown inner state variable and hardly accessible 
by measurements. Hence, an observer is employed to estimate 
this pressure using the measured pressure at the outlet throttle of 
the fuel cell system. Thereby, the dynamic observer allows for 
detecting temporary shortages or peak pressures at the electrode 
and enables the control system to counteract accordingly. 


2. Fuel cell subsystem 


The core of the test rig is given by a fuel cell stack as depicted 
in Fig. 3 from the centre of solar energy and hydrogen research 
(ZSW) in Ulm with 25 cells and a nominal power of 800 We1. The 
air supply subsystem of the fuel cell comprises the flow chan- 
nels of the fuel cells and all the peripheral components required 
to supply air and to remove water and exhaust. The test rig is 
supplied with dry compressed air with a maximum of 2 baraps. 
The air is humidified by a dew point humidifier. The air flow fed 
to the stack is controlled by a mass flow controller. The supply 
manifold between the MFC and the stack consists of pipes and 
connections. Another manifold connects the outlet of the stack 
with the outlet throttle that is used for controlling the cathode 
pressure. 


Fig. 3. Fuel cell stack. 
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Fig. 4. Structure of the gas supply subsystem. 


In the following, the supply pressure at the MFC is assumed 
as constant. Modelling focuses on the system illustrated in the 
block diagram in Fig. 4. 


3. Dynamic model of the air supply system 


Modelling the air supply system aims at describing the inter- 
action of the mass flow controller and the outlet throttle with the 
cathode pressure pc and the excess ratio of oxygen ào, dynami- 
cally. The model is achieved by developing an equivalent circuit 
diagram for the pneumatic system using analogies to electrical 
systems [10]. 


3.1. Lumped pneumatic model 


Pneumatic resistances and capacities form lumped parame- 
ters of the model. A pneumatic resistance describes the friction 
losses of a gas flow and a capacity represents the ability of 
a component to store mass [9]. Pipes, manifolds, throttles 
and flow channels of the fuel cell stack can be modelled by 
using such pneumatic resistances and capacities. In the follow- 
ing these subsystems as well as the whole equivalent circuit 
including the mass flow controller will be described mathemat- 
ically. 


3.2. Model of the mass flow controller 


As described in (1) the desired air mass flow to operate the 
system is calculated using Faraday’s Law 


ào, N Mo, [rc 


(4) 
yo, z F 


MMFC,d = 
Since the MFC does not react instantaneously on a change in the 
desired mass flow 7impc.a, this characteristic can be described 
as a first order lag element with the time constant t 


d 
T MFC + MFC = MMFC,d. (5) 


If the pressure pı at the outlet of the MFC rises, the effective 
mass flow mı decreases. If pı equals the supply pressure pmax 
the mass flow saturates which is described by the saturation 
function 


m; = tanh{a (pı — Pmax)}/™MFc- (6) 
3.3. Model of the manifold 


There are two manifolds in the system. The inlet manifold 
between the MFC and the fuel cell and the outlet manifold 


Min Mout 


Pin Pout 


Fig. 5. Model of the manifold. 


between the fuel cell and the throttle. The simplest model of 
a manifold or a pipe is a pneumatic capacity and a pneumatic 
resistance in parallel (Fig. 5). Generally, there are two borderline 
cases of flow, laminar [10] and turbulent [4] flow. Accordingly, 
laminar resistances Rj, and turbulent resistances Riyp have to 
be distinguished. In the case of laminar flow the mass flow is 
proportional to the pressure drop at the resistance 


Ap = Ramm. (7) 


In the case of turbulent flow the mass flow is proportional to the 
root of the pressure drop 


V Ap = Kub, (8) 


with the constant parameter Kurb and the pressure dependent 
turbulent flow resistance 


Rturb = Kturb y Ap. (9) 


Since laminar and turbulent flow are borderline cases, the 
pneumatic resistance of a pipe can also be modelled as a series 
of a laminar and a turbulent resistance (Fig. 6). 


3.4. Model of the fuel cell stack 


In analogy to the manifold model the fuel cell stack model 
(Fig. 7) consists of a laminar and a turbulent resistance and a 
capacity. The resistance results from the friction losses of the 
gas flow through the flow fields. The capacity results from the 
storage volume of the flow fields and the gas diffusion layer. 
The FC model differs in the configuration of the components, 
which is founded on the results of the parameter identification 
(cf. Section 4). Furthermore, the change in mass flow at the 
electrode has to be taken into consideration. The effective mass 
flow 


Metf = My + Mmem + MHO (10) 


Pin | Pout 


Fig. 6. Equivalent circuit of the inlet and outlet manifold. 
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Fig. 7. Equivalent circuit of the fuel cell stack. 


is modelled as a source in the equivalent circuit. The effective 
mass flow consists of the mass flow m, of the reaction at the 
cathode side, the mass transport mem through the membrane 
and the mass flow 74,0 of water coming from the humidifier. 
The mass flow of the reaction 


jee ee wo: NIgc 
i ZH20 20> F 


(1) 


comprises the creation of water and the consumption of oxygen 
at the current Ipc. The mass flow 

Mmem = Mdrag — Mdiff (12) 
describes the water transport through the membrane caused by 
the electro-osmotic drag (7itgrag) and by back-diffusion (maite). 
The description of marag and maige can be found in Refs. [7] 
and [2]. The derivation of the amount of water mmy,0 car- 
ried in by the humidifier can be seen in Ref. [13]. It results 
in 

Rair Psat(Thum) 


myo = (13) 
i Ry,0 PC — Psat(Thum) 


with the gas constants Rair, Rmo and the saturation pressure 
Psat at the temperature Thum of the humidifier. 


3.5. Model of the throttle 


The throttle at the outlet is modelled by a turbulent resistance 
[4]. The throttle resistance Rr can be changed by varying the 
opening area Ar. The mass flow through the throttle is calculated 
by 


Tow = ee = pA any? Jpeg. (14) 


3.6. Model of the air supply system 


Connecting all subsystems to the overall model of the air 
supply system leads to the equivalent circuit in Fig. 8 with the 
combined resistances 


Ry = Riami + Rturbi + Riam2; 


Ry = Rub2 and Rm = Riam + Rturb3- 


my ms m Mout 
o> > + > + 
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Fig. 8. Equivalent circuit of the air supply system. 


The dynamics of the pneumatic system is described by a set of 
oridnary first order differential equations 


Pı — P2 


mı =C ppt Ri (15) 
1— p2 , 2—- v 
ay a age (16) 
Ry 
p2 — p3 , p3 
—— = C3 p3 + ——_. (17) 
Ry á Rig + Rr 


These differential Eqs. (15)—(17) together with the equations of 
the MFC (5) and (6) describe the entire system dynamics of the 
air supply system. 


4. Parameter identification 


In Ref. [13] measurements were done to identify the parame- 
ters C1, C2, C3, Rī, Ry and Ry of the air supply model for a wide 
range of mass flows. The air mass flow through the system was 
changed stepwise from a minimum of 0.5 gs7! to a maximum 
of 2gs~!. The parameters for each subsystem were obtained 
from the measurements by using the Nelder—Mead algorithm 
for parameter optimisation [5]. Inputs to the algorithm are the 
measured pressure values Yin , the actual mass flow min and 
the derived model description i in form of discretised differen- 
tial equations. With these finite difference equations the model 
values Vinod 3E calculated. The algorithm then calculates the 
optimal parameter vector 


T 
0 = [C1, C2, C3, Riami, Riam2, Riam3, Kturb1, Kturb2, Kturb3] 


by minimizing the quadratic error 


C= Oia Vmod) (18) 


between the measured values and the model values. 

As an example of an identified subsystem, in Fig. 9 the results 
from the parameter identification for the fuel cell stack can be 
seen. Here the measured pressure values p2 mea and the simu- 
lated model values p2 moa are compared. 

By parameter identification the values of the parameters for 
the MFC (rt), the capacities (C1, C2, C3), the laminar resistances 
(Riam1, Riam2, Riam3) and the constant factors for the turbulent 
resistances (Kturb1, Kturb2, Kturb3) were determined as can be 
seen in Table 1. 

The last missing parameter in Fig. 8 is the throttle resistance 
Rr. As described above the throttle is modelled by a turbulent 


M.A. Danzer et al. / Journal of Power Sources 176 (2008) 515-522 519 


0 50 100 150 
t/s 


Fig. 9. Comparison between measured values for p2 and model values. 


resistance. By combining (8) and (9) you get 
Rr = K7yritow (19) 


The factor Kr depends on the throttle opening area Ay, which 
can be changed by varying the control voltage Vr of the throttle. 
The throttle resistance in (19) depends on the mass flow Mout 
through the throttle and the control voltage Vr. In the model 
the value for the throttle resistance Rr is derived from a char- 
acteristic map Rr = f(Mou, Vr), which can be seen in Ref. 
[13]. 


5. Nonlinear model-based control 


The main control objectives are obtaining high overall effi- 
ciency and avoiding starvation. As for the air supply system of 
the fuel cell, the multivariable control involves set-point con- 
trol of the cathode pressure pc = p2 and set-point control of 
the oxygen excess ratio A49,. Accordingly, the control outputs 
are chosen as yj = p2 and y2 = Ag,. The subsystem offers two 
means that have an effect on the pressure and the excess ratio: the 
MFC at the inlet and the outlet throttle. Thus, the control inputs 
are the desired mass flow rca and the flow resistance Rr of 
the throttle, which refers to the cross sectional area (cf. (14)) 
and the control voltage, respectively. Here, the integrated MFC 
represents a fast underlying control loop. The state variables of 


Table 1 

Identified parameters 

Parameter Value 

T 0.495 s 

Ci 2.17 x 107°? ms? 
C2 4.46 x 1078 ms? 
C3 1.25 x 107° ms? 
Rigi 6.94 x 10° (ms)! 
Riam2 1.53 x 107 (ms)! 
Riam3 1.43 x 10° (ms)! 
Kurbi 47,940 (mkg)~°> 
Kiurb2 54,901 (mkg)~°° 
Kiurb3 44,651 (mkg)~°5 


stabilising | V inverse | U flat y 
input dynamics system 


4 


Fig. 10. Structure of the flatness-based control. 


the system are given by the pressures at the capacities: 


; pı 
MMFC,d = P2 
u= >, x=]|]p|, Y= . (20) 
Rr ~ ho» 
P3 


As the system is nonlinear, flatness-based control techniques are 
envisaged. 


5.1. Principles of flatness-based control 


Differential flatness is a geometrically motivated concept for 
analysis and synthesis of nonlinear systems. According to Ref. 
[8], nonlinear flat systems are equivalent, in a more general way, 
to linear controllable systems (Fig. 10). 

A system is denoted flat if an output vector y = [y1, ... 
exists that fulfills the following conditions: 


3 Yp 


(i) The output variables y; can be stated as functions of the 
state variables x;, the input variables u; and afinite number 
of their time derivatives 


y = D, u, ù, ..., u). (21) 


(ii) All state variables x; and all inputs u; can be stated as 
functions of the output variables y; and a finite number of 
their time derivatives 


x= DQ, j, 9), (22) 


u = (y, Peoceg yh): (23) 
(iii) The dimension of the input vector is equal to the dimension 
of the output vector 


dimu = dimy. (24) 


If the output variables of interest can be proven flat outputs, 
trajectory tracking control design becomes straight-forward. The 
dynamics of the resulting linear error dynamics can be specified 
by introducing a new stabilising input 


(6+1) 


vi = year? + kph yD +--+ + koa — yi) (25) 


and replacing the highest derivative of y; in the expression of 
the control input with v; according to (25). This results in the 
following inverse dynamics 


u = Waly, ¥,-.., yf, V), (26) 


where the control inputs are calculated depending on the mea- 
sured and desired outputs. 
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5.2. Application to the system 


For the proof of flatness the system description from Section 
3.6 is examined. Unfortunately, the proof of the flat output can- 
didate ý according to (20) fails. Instead, it can be shown that 
y= [yı ya]? = [pı um]! represents a flat output. The desired 
value pı a for the inlet pressure can be derived from the steady- 
state values of the desired value p2, for the cathode pressure, 
the actual mass flow mı and the flow resistance Ry 


Pid = prat Ry. (27) 


The actual values of the turbulent resistances are calculated by 
using the steady-state values of the according pressure differ- 
ences. The relative degree of the first input yı is three, whereas 
for the second input y2 the relative degree becomes one. 
According to (25) and (26), the new stabilising input 


TE Yı | _ | Yra + kizëêi + kii éi + 10 e1 (28) 
y2 yaa + k20 e2 

is introduced and inserted into the inverse dynamics 
m V21 (Y1, V1, Ÿ1, Vi, Y2, V 

u= | MBCA | _ a1 V1, Y1, V1, Y2, v2) l (29) 
Rr W0(y1, Y2, V2) 

Thus, the linear tracking error dynamics 

pi k .. k . k = 

e1 +k +k ei + kine: = 0 (30) 


ér + kx e2 = 0,7 


with ei = yia — yi can be specified by the coefficients k12, k11, 
kio and ky. 


6. Cathode pressure observation 


The aim of the feedback control is to maintain a constant cath- 
ode pressure. However, the pressure pz at the electrode inside 
the fuel cell stack is a state variable that cannot be measured 
directly. The only measurable state variable is the pressure p3 at 
the outlet throttle. Since the turbulent flow resistances increase 
with raising pressure, it is not advisable to control the output 
pressure p3 instead of the cathode pressure p2. Especially at 
high mass flows both pressures differ significantly (cf. Fig. 16). 
An alternative is to estimate the cathode pressure dynamically. 

The fuel cell current is measured and therefore, the actual 
utilisation follows from Faraday’s Law (1). The actual value of 
the inlet mass flow mı is available from the MFC. The actual 
flow resistance Rr of the throttle is known as well. Under these 
preconditions the pressure p2 at the electrode can be estimated 
by a tracking observer. 

For the observer design the input u, the measured output ym 
and the states x are chosen as follows 
Ym = P3 (31) 


u=tn, x=(pi, p2, pal’, 


and lead to the corresponding nonlinear state space description 


k= f(x, u), (32) 


flat 
system 


stabilising | V inverse 
input dynamics 


system 
model 


Fig. 11. Structure of the nonlinear tracking observer. 


y = A(x) = x3. a 


Consequently, the flatness-based feedback control (29) can be 
realised with the estimated state vector * provided by a nonlinear 
tracking observer 

x= fQ, @) +A Om — A), (34) 


with time-varying observer-gain L(t) [8] (Fig. 11). 

The derivation of the observer-gain is discussed in detail in 
Ref. [3]. For the air supply system the time-varying observer- 
gain results in 


Ut) = AA, kO, BOM, (35) 


where all components depend on the pressure dependent tur- 
bulent resistances and the third gain additionally on the actual 
throttle resistance Rr. 


7. Simulation results 


In the following, the effectiveness of the proposed nonlinear 
multivariable control combined with the observer is investigated 
by simulations, where the fuel cell current acts as disturbance on 
the fuel cell system. The current variation is given by a sequence 
consisting of a steep increase from 0.1 A cm7*to 0.4Acm7?, 
a subsequent constant level of 0.4 A cm7?, and a steep decline 
to 0.1 A cm7? as shown in Fig. 12. 

The resulting time response of the control outputs, i.e., the 
excess ratio of oxygen Ag, and the cathode pressure p2, are 
depicted in Figs. 13 and 14, respectively. 

During steady-state, both excess ratio of oxygen Ao, and the 
cathode pressure p2 match their desired values. At the rising 
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Fig. 12. Fuel cell current. 
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Fig. 13. Excess ratio of oxygen Ao, using observer-based multivariable control. 


edge, the excess ratio reaches a minimum of 1.9 and a maxi- 
mum of 4.6 at the falling edge. Since the excess ratio never falls 
below 1, the fuel cell is never in danger of starvation. Due to the 
increase in current it is expected that the pressure level at the 
cathode decreases. But as can be seen in Fig. 14, p2 increases 
slightly from 5 x 10+ to 5.15 x 10* Pa. This has to be seen 
in the context that pı is controlled with the input variables in 
(29). Due to diffusion effects through the gas diffusion layer the 
occurring pressure increase in the cathode is advantageous at 
rising currents. 

The observed values of the cathode pressure p2 follow 
the simulated values with high accuracy. The corresponding 
observer error €o,p, = P2 — p2 of the anode pressure is depicted 
in Fig. 15. 

To emphasise the effectiveness of the proposed flatness-based 
control extended with the nonlinear tracking observer, a com- 
parison with the commonly used combination of feedforward 
control of the mass flow rc and PI-control of the outlet pres- 
sure Pout at the outlet throttle has been performed. Here, the same 
current trajectory as depicted in Fig. 12 has been employed. 

Since the outlet pressure Pout is controlled instead of p2, the 
cathode pressure is subject to steady-state error, especially at 
high current densities (cf. Fig. 16). At the first load change the 
cathode pressure drops and at the second load change to lower 
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Fig. 14. Cathode pressure pz using observer-based multivariable control. 
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Fig. 15. Observer error of the cathode pressure. 
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Fig. 16. Cathode pressure p2 with the commonly used PI-control of the outlet 
pressure Pout- 


current values a peak pressure is observed inside the fuel cell 
stack. 

The feedforward controlled excess ratio of oxygen reaches 
values close to 1 at the rising edge of the current trajectory (cf. 
Fig. 17). For this period, the fuel cell runs the risk of starva- 
tion. 


0 5 10 15 20 


Fig. 17. Excess ratio of oxygen Ào, with the commonly used feedforward 
control of the mass flow. 
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8. Conclusion and outlook 


As basis for both control and observer design, a dynamic 
model of the air supply subsystem at the cathode side of the fuel 
cell is derived. The control objectives in terms of obtaining high 
overall efficiency and avoiding starvation are achieved by mul- 
tivariable control of both cathode pressure and excess ratio of 
oxygen. With the designed control structure the fuel cell system 
can be operated dynamically at lower excess ratios of oxy- 
gen without running the risk of starvation. The non-measurable 
anode pressure is estimated by a nonlinear tracking observer. 
The dynamic pressure observer allows for detecting pressure 
drops and peak pressures in the cathode volume and, thereby, 
reduces the danger of starvation and overpressure, respectively. 

The effectiveness of the control system is shown by com- 
paring the simulation results of the proposed nonlinear control 
structure with the results of a usually used feedforward control of 
the mass flow combined with a PI-control of the outlet pressure. 
At this, the same current trajectory is employed. 

Future work envisages a system where the MFC is replaced 
by a compressor. In the model the dynamics of the MFC is 
substituted by a model of the compressor dynamics. Finally, 
an air supply system is obtained with two voltages as control 
variables, one for the compressor and the second for the outlet 
throttle. To realise the same multivariable control and observer 
design as discussed above two pressure sensors are sufficient: 
one at the outlet of the compressor to determine the mass flow 
into the system and the second, as hitherto, at the outlet throttle. 

To further minimise the danger of starvation the developed 
control system can be extended by aload governor [11]. Thereby, 
the current drawn from the fuel cell can be constricted depending 
on the observed pressure at the electrodes, so that the constraint 
on the oxygen partial pressure inside the cathode is strictly 
enforced. 
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